Effects of partial thermalization on HBT interferometry 
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Abstract 

Hydrodynamical models have generally failed to describe interferometry radii measured at RHIC. 
In order to investigate this "HBT puzzle", we carry out a systematic study of HBT radii in ul- 
trarelativistic heavy-ion collisions within a two-dimensional transport model. We compute the 
transverse radii R a and R s as functions of p, for various values of the Knudsen number, which 
measures the degree of thermalization in the system. For realistic values of the Knudsen number 
estimated from V2 data, we obtain R /R s - 1.2, much closer to data than standard hydrodynam- 
ical results. Femtoscopic observables vary little with the degree of thermalization. Azimuthal 
oscillations of the radii in non central collisions do not provide a good probe of thermalization. 



1. A short reminder on the model 

In this talk, we review our recent calculation [1] of HBT radii, using a transport model de- 
scribed in more details in 12J] . We model the expansion of the matter created in the central rapidity 
region of a heavy-ion collision using a relativistic 2 dimensional Boltzmann equation with mass- 
less particles. The equation of state (EOS) of the medium is that of a conformal ideal gas s = 2P. 
We use two different initial momentum distributions. The first is a locally thermal distribution 
in order to make comparison with hydrodynamic simulations [2]. The second is a more realistic 
distribution based on a parametrisation of the Color Glass Condensate (CGC) initial gluon spec- 
trum from Q3l |4£| . In both cases, our momentum distribution is scaled such a way that the average 
transverse momentum per particle is given by (p,) = 420MeV, which corresponds roughly to the 
value for pions at the top RHIC energy. 

In our numerical solution of the 2+1 dimensional Boltzmann equation, the partonic cross 
section is assumed to be constant and isotropic in the center of mass frame for simplicity. The 
dynamics of the system is then controlled by a single quantity, the mean free path A, which is the 
average distance travelled by a particle between two collisions. The Knudsen number K is the 
ratio of the mean free path to the characteristic size of the system. Because K^ 1 is proportional 
to the average number of collisions per particle during the evolution, we use K to characterize 
the degree of thermalization of the system; hydrodynamics is the limit K — > while the limit 
K — » oo corresponds to the free-streaming regime. 

For a particle with outgoing momentum p ( , we denote by (f, x,y) the space-time point where 
the last collision occurs. The "out" and "side" coordinates are then defined as the projections 
parallel and orthogonal to the particle momentum. The HBT radii are then obtained from rms 
widths of the space time distribution of particles with the same outgoing momentum. Radii 
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defined in this way coincide with those obtained from the curvature of the correlation function at 
zero relative momentum. 

2. Central collisions 

In order to mimic a central Au-Au collision at RHIC, we use a gaussian initial density profile 
with cr x — cr v — 3 fm. This value corresponds to the rms width of the initial density profile in an 
optical Glauber calculation. 
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Figure 1: HBT radius R a versus transverse momentum p, of particles: (left) Transport calculations; the curves are labeled 
by the value of the Knudsen number K. (right) Ideal hydrodynamics at fixed freeze out time t; the curves are labeled by 
the value of t. 

Figure[T](left) displays R a versus p, for different values of the Knudsen number. It shows that 
R a decreases when increasing K 1 or p,. The decrease with p, is usually considered as a signature 
of the hydrodynamical behaviour, but even if it is more pronounced for hydrodynamics, it is also 
seen in the free streaming regime. Figure [T] (right) displays R„ versus p, for different values of 
the freeze-out time in an ideal hydrodynamic calculation. Comparing figure Q] (left) and (right), 
we see that hydrodynamics at t — 0.2fm/c gives the same radii as transport for high values of K 
and that hydrodynamics at t — 40fm/c gives almost the same results as transport calculations for 
small K. Our interpretation is that decreasing K amounts to increasing the number of collisions 
per particles until the "freeze out", which has qualitatively the same effect as increasing the 
duration of the hydrodynamic phase. 

In the transport simulation, we observe that R„ increases and R s decreases when the number 
of collisions per particle l/K increases, as plotted on figure [2] (left). On the same figure, we also 
plot the elliptic flow from J2[] versus K 1 . One sees that V2 saturates to the hydro limit much 
faster than R„ and R s . Hydrodynamical calculations Jit] |0] J3l usually yield a value of R a /R s of 
the order of 1.5, while RHIC data are compatible with 1. Figure [2] (right) displays the HBT ratio 
R /R s with K = 0.3, which is the value obtained by fitting the centrality dependence of elliptic 
flow |g]. Our results are in much better agreement with experimental data than models based 
on ideal hydrodynamics. It has been recently argued Jioll that ideal hydrodynamics simulations 
with an early freeze out also explains the HBT puzzle. Our results also show that increasing the 
Knudsen number amounts to decreasing the freeze out time. We find that partial thermalization, 
which has been shown to explain the centrality dependence of V2, also solves most of the HBT 
puzzle. In an other recent work, Scott Pratt argued that the HBT puzzle results from a 
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Figure 2: Results obtained for HBT radii averaged over the p, interval [0.25 ,0.75] GeV: (left) R„ and R s versus the 
degree of fhermalization. The dotted curve shows, for sake of illustration, the variation of elliptic flow in a non-central 
collision, scaled by the hydrodynamical limit (from |2]). (right) R /R s versus p t for K = 0.3. Dashed line: thermal initial 
conditions; full line: CGC initial conditions. Stars: data from STAR |9] 



combination of several factors. Our results show that partial thermalization alone solves most of 
the puzzle. 



3. Azimuthally sensitive HBT 

For a non-central collision, the interaction region is elliptic, and HBT radii depend on the 
azimuthal angle </>, defined with respect to the impact parameter direction. Figure [3] (left) shows 
the radii in plane and out of plane as functions of the average number of collisions per particle. As 
in the case of central collisions, we observe that R increases and R s decreases when increasing 
l/K, but the slope is larger for <p — than for <p = n/2. Our interpretation is that thermalization 
is faster in plane that out of plane due to flow effects. 
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Figure 3: Evolution of AzHBT observables with the degree of thermalization: (left) In-plane (ip = 0) and out-of-plane 
((j> = tt/2) radii for thermal initial conditions. (right) Ratios of oscillation amplitudes for thermal initial conditions. 



We define three oscillations amplitudes: 

ARl = R 2 o (k/2)-R 2 o (0) 
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AR 2 = R 2 s (0)-R 2 s (n/2) 
AR os = R os (3n/4)-R os (n/4). 



(1) 



These amplitudes scale with the eccentricity of the overlap area, which is not directly measured. 
We avoid this dependence by taking ratios as shown on figure [3] (right). These ratios are equal to 
unity in the free streaming regime and can also be measured experimentally. For realistic values 
of K in the range [0.3, 0.5], both ratios deviate little from unity. On figure [3] (right) we also plot 
the ratio of the eccentricity seen in HBT radii e s = (R 2 (0) - R 2 (tt/2))/(R 2 (0) + R 2 (tt/2)) to the 
initial eccentricity. This ratio also remains close to unity. We thus conclude that none of these 
observables may be a good probe of thermalization. 

Although our model is too crude to reproduce the absolute magnitude of HBT radii, we 
expect that the above ratios are less model dependent. Our results for AR os /ARl and tsR 2 ltsR 2 
are compatible with experimental data [1], but the latter have large error bars. On the other 
hand, we miss the value of e s , which seems to be washed out by the expansion. Hydrodynamical 
calculations using a soft EOS, reported a fair agreement with the mesured value of e s . It thus 
seems that the soft equation of state of QCD is responsible for the reduced eccentricity seen in 
the data. 

Conclusion 

We have carried out a systematic study of how HBT observables evolve with the degree of 
thermalization {K). Our results show that HBT observables depend only weakly on K. The 
decrease of the radii with p, is already expected from initial conditions and is only slightly 
enhanced by collective flow. We have also shown that for realistic values of K, inferred from 
elliptic flow study, the ratio R /R s is lower than 1.2. Partial thermalization solves most of the 
HBT puzzle. 
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